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Abstract 

We study a Brownian motor driven by a colored non Gaussian noise source with a q- 
dependent probability distribution, where q is a parameter indicating the departure from 
Gaussianity. For q — 1 the noise is Gaussian (Ornstein-Uhlenbeck), while, for q > 1, 
the probability distribution falls like a power law. In the latter case, we find a marked 
enhancement of both the current and the efficiency of the Brownian motor in the overdamped 
regime. We also analyze the case with inertia and show that, again for q > 1, a remarkable 
increase of the ratchet's mass separation capability is obtained. 
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The study of noise induced transport by "ratchets" has attracted an increasing number 
of researchers that have produced a vast literature [[j], 0, ||, [|, |5|. Among other aspects, the 
motivation of these studies has been prompted by both their possible biological interest Q as 
well as their potential technological applications. Since the initial works, besides the built-in 
ratchet-like bias and correlated fluctuations (see for instance Q), different aspects have been 
studied like tilting @ and pulsating potentials, velocity inversions ||, etc. There are some 
reviews where it is possible to grasp the state of the art ]|, . 

So far, almost all studies have used Gaussian noises, with few exceptions that mainly ex- 
ploited dichotomic processes @ . Here we analyze the effect of a particular class of colored non 
Gaussian noise on the transport properties of Brownian motors. Such a noise source is based on 
the so called Tsallis statistics || with a probability distribution that depends on q, a parameter 
indicating the departure from the Gaussian (q = 1) behavior. Some of the motivations for study- 
ing the effect of non Gaussian noises are, in addition to the intrinsic interest within the realm 
of noise induced phenomena, different experimental data indicating that for some biologically 
motivated systems fluctuations can have a non Gaussian character. An example are current 
measurements through voltage-sensitive ion channels in a cell membrane or experiments on the 
sensory system of rat skin Q. It is also worth remarking here that recent detailed studies on 
the source of fluctuations in some biological systems [[h]] clearly indicate that noise sources in 
general could be non Gaussian and their distribution bounded. Here, again as a consequence 
of the non Gaussian character of the driven noise, we find a remarkable increase of the current 
together with an enhancement of the motor efficiency that, additionally, shows an optimum for a 
given degree of non Gaussianity. Moreover, when the inertia is taken into account we find that, 
when departing from the Gaussian case, there is a remarkable increment in the mass separation 
efficiency of these devices. 

We start considering the general system 

™§ = "7§ - V'{x) -F + ^(t) + v(t), (1) 

where m is the mass of the particle, 7 the friction constant, V(x) the ratchet potential, F 
is a constant "load" force, and the thermal noise satisfying (£(£)£(£')) = ^KT5{t — t'). 
Finally, rj(t) is the time correlated forcing (with zero mean) that keeps the system out of thermal 
equilibrium allowing the rectification of the motion. For this type of ratchet model several 
different kinds of time correlated forcing have been considered in the literature || 0. We 
may distinguish the cases of deterministic (such as the periodic forcing and the "synthetic 
noise" discussed in the pioneering work of Magnasco ||) and of stochastic forcing. Among the 
stochastic forcings the more usual noise sources are dichotomic and Ornstein-Uhlenbeck (OU) 
noises ]?]]. Here, as indicated above, we will study the effects of a different family of stochastic 
forcing having non Gaussian statistics. The main characteristic introduced by the non Gaussian 
form of the forcing we consider here is the appearance of arbitrary strong "kicks" with relatively 
high probability when compared, for example, with the OU Gaussian process. As we shall see, 
in a general situation (without fine tuning of the parameters) , this leads to the above indicated 
remarkable effects. 

We will consider the dynamics of rj{t) as described by the Langevin equation [[Tl|, 12] 



with <C(t)C(*)'> = 2ZM(i ~ ?) and 



".W-^ijMi + sh-D^]. (3) 
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Previous studies of such processes in connection with stochastic resonance problems |TT| , [l2| and 
dynamical trapping [^] have shown that the non Gaussianity of the noise leads to interesting 
effects. For q = 1, the process rj is the (Gaussian) Ornstein-Uhlenbeck one (with correlation time 
equal to r), while for q ^ 1 it is a non Gaussian process. For q < 1 the probability distribution 
has a cut-off at u> = [(1 — q)r/(2D)]~2 , and is given by: 




P q <l(v) = \ ^ 'f " (4) 

otherwise 

where Z q is the normalization constant. For 1 < q < 3, the probability distribution (for — oo < 
rj < oo) is 



P q >i(v) = ^ 



1 + 



1 



(5) 



While keeping D constant, the width or dispersion of the distribution increases with q. This 
means that, the higher the q, the stronger the "kicks" that the particle will receive. The second 
moment of the distribution (or intensity of the non Gaussian noise, that we call D ng ) diverges 
for q > 5/3 ~ 1.66 while, for q < 5/3, is given by 

Here, and in what follows, we shall keep q < 5/3. 

For the ratchet potential we will first consider the same form as in [|j (with period 2tt) 

V(x) = — j dx(exp[a cos(x)]/ Jo(ia) — 1), (7) 

with a = 16. The integrand is the ratchet force (-V'(x)) appearing in Eq.(|l|). 

Firstly, we will analyze the over damped regime setting m = and 7=1. We are interested 
on analyzing the dependence of the mean current J = (^|) and the efficiency e on the different 
parameters. In particular, their dependence on q, the parameter indicating the degree of non 
Gaussianity of the noise distribution. The efficiency of the ratchet system is defined as the ratio 
of the work (per unit time) done by the particle "against" the load force F, into the mean power 



injected to the system through the external forcing 7/ [14] 



For the numerator we get F(^-) = FJ, while, for the denominator 

Interesting and complete discussions on the thermodynamics and energetics of ratchet systems 
can be found in fj"5]| . 

It is worth mentioning here that the "effective Markovian approximation" introduced in [|l^] 
is not adequate for the present case as it only properly works in the low r limit. Here, in the 
overdamped regime we are able to give an approximate analytical solution for the problem, 
which is expected to be valid in the large correlation time regime (-^ >> 1): we perform the 
adiabatic approximation of solving the Fokker-Planck equation associated to Eq. (|l|) assuming 



a constant value of r/ [16]. This leads us to obtain a value for the current J(r/) which depends 
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on i] and, in order to get the final result, we should perform an average of J(rj) over rj using the 
distribution P q (r)) with the desired value of q. 

In Fig. 1, we show typical analytical results for the current and the efficiency as functions 
of q together with results coming from numerical simulations (for the complete system given 
by Eqs.(jl]) and (||)). Calculations have been done in a region of parameters similar to the one 
studied in Q but considering (apart from the difference provided by the non Gaussian noise) 
a non-zero load force that leads to a non-vanishing efficiency. As can be seen, although there 
is not a quantitative agreement between theory and simulations, the adiabatic approximation 
predicts qualitatively very well the behavior of J (and e) as q is varied. As shown in the figure, 
the current grows monotonously with q (at least for q < 5/3) while there is an optimal value 
of q (> 1) which gives the maximum efficiency. This fact is interpreted as follows: when q is 
increased, the width of the P q {rj) distribution grows and high values of the non Gaussian noise 
become more frequent, this leads to an improvement of the current. Although the mean value of 
J increases monotonously with q, the grow of the width of P q (rj) leads to an enhancement of the 
fluctuations around this mean value. This is the origin of the efficiency's decay that occurs for 
high values of q: in this region, in spite of having a large (positive) mean value of the current, 
for a given realization of the process, the transport of the particle towards the desired direction 
is far from being assured. 

In Fig. 2 we show results from simulations for J and e as functions of q for different values 
of D, the intensity of the white noise in Eq. (|2|). The results correspond to KT = 0, hence, 
the only noise present in the system is the non Gaussian one. On the curve corresponding to 
the results for J (Fig. 2. a.), we indicate with error bars the dispersion of the results. The huge 
growth of the dispersion occurring for q > 1.3, as well as the decay of the efficiency for the same 
values of q, is apparent. 

It is worth recalling here that q = 1 corresponds to the Gaussian OU noise (analyzed, for 
example in [Q] and [17]]). Hence, our results show that the transport mechanism becomes more 
efficient when the stochastic forcing has a non Gaussian distribution with q > 1. 

In the previous calculations, we have analyzed the values of J and e as functions of q for 
fixed values of D. However, there are situations where we should consider that the non Gaussian 
noise source is the "primary" source, for instance see [10] for biologically motivated problems. 
Besides this, there could also be situations of clear technological interest. For these reasons we 
wondered which would be the dependence of J and e on q when we fix the intensity of the non 
Gaussian noise D ng defined in Eq.(^). In such case, the energy per unit time supplied to the 
ratchet is independent of q. Our results (not shown) indicate that, at variance to what occurs in 
calculations for fixed D, for constant D ng , the optimum value of q which maximizes the efficiency 
also maximizes the mean current J, which decays for larger values of q. This result is easy to 
understand by observing that the width of the P q {ri) distribution decreases with q when D ng is 
kept constant. 

Now we turn to study the m ^ case, that is, the situations in which the inertia effects 
are relevant. In Fig. 3 we show the dependence of the current J on the mass m for different 
values of q. The results are from simulations for zero temperature and without load force. It 
can be seen that, as m is increased from 0, the inertial effects initially contribute to increase 
the current, until an optimal value of m is reached. As it is expected, for high values of m, the 
motion of the particle becomes difficult and, for m — > oo, the current vanishes. 

An interesting effect appears for q = 1.3: for a well defined interval of the value of the mass 
(ranging approximately from m = 100 to m = 5000), a negative current is observed. This can 
be explained as a consequence of the high value of the mass, that makes the inertial effects 
much more important than those of the ratchet potential. Due to the high value of m, most 
of the time the particle will be almost at rest around the minimum of the ratchet potential. 
However, a very high value of r/ may occur (that may last for an interval of time of the order 
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of r) causing the particle to move a period of the potential to the right or to the left. The 
obvious question is, what is the value of the force r] needed in order to make the particle jump 
a period of the potential to the left (right), during a time r?. Assuming no friction (which is 
of no importance for high values of m), a constant (average) value of the ratchet force V , and 
performing a classical calculation we get: 

, /, 2md , . 

1= V +-2-, 10 

where d is the distance from one minimum of V{x) to the next left (right) maximum of V(x). In 
our system, for a jump to the left, we have \ V'\ ~ 4.5 and d ~ 1.1, while, for a jump to the right, 
we have \V'\ ~ .96 and d ~ 5.2. Considering r = 100/(27r), for m = 800 (the value for which 
in Fig. 3 the minimum of the current is obtained for q = 1.3) we find that the required value 
of r\ for a jump to the left (right) is \rji\ ~ 6.9 (\rj r \ ~ 32.4). Hence, for this value of the mass, 
the jumps to the left are more probable than the jumps to the right, and a negative current is 
to be expected. In contrast, for a value of the mass m = 100 (approximately where the current 
change from being positive to being negative) we find \r]i\ ~ 5.4 and |ry r | ~ 5.1, i.e. both jumps 
are almost equiprobable. Note that, for m = 800, the large value of m makes the second term 
in Eq. ( |To[ ) more important than the contribution of \V'\ (the inertial effect dominates over the 
ratchet force), hence, the difference in the value of r] necessary for a jump to the left or to the 
right comes, essentially, from the difference on the distances (d) that the heavy particle should 
be displaced to each side. 

Until this point, the analysis makes no mention of the value of q of the noise distribution. The 
fact that the (J < 0)~effect appears for q = 1.3 and not for q = 1 is understood as a consequence 
of the fact that, for q = 1, the occurrence of a value of noise \tjl\ is highly improbable in the 
regime where the inertia dominates over the effect of the ratchet force. Hence, the particle 
remains essentially motionless. 

These results imply that separation of masses (particles with different masses moving in 
opposite directions) occurs, and that this happens in the absence of load force, due to the "big 
help" given by the non Gaussian noise. To our knowledge, separation of masses by ratchets 
with zero load force have not being reported before. However, the separation found here occurs 
for particles with a ratio of masses of the order of 10 or more (say 800/80 in Fig. 3.), while it 



has been shown in 18, 19| that, with a load force and considering simply OU noise, particles of 
much closer masses can be separated by ratchets. 

Now we analyze in more detail the problem of mass separation, an aspect that has been 
studied in some works ||, U|. In 

view of the results discussed above, it is reasonable to 
expect that non Gaussian noises may improve the capability of mass separation in ratchets in 
more general situations. Reference [O] was one of the primary works discussing mass separation 
by ratchets. There, the authors analyzed a ratchet system like the one described by Eq. (|l]) 
considering OU noise as external forcing (in our case it corresponds to q = 1). They studied 
(both numerically and analytically) the dynamics for different values of the correlation time of 
the forcing r, finding that there is a region of parameters where mass separation occurs. This 
means that the direction of the current is found to be mass-dependent: the "heavy" species 
moves in the negative sense while the "light" one, do so in the positive sense. 

Here, in order to compare results, we analyze the same system studied in [18, 19| but 
considering the non Gaussian forcing described by Eq.(|l]). Hence, we study the system in 
Eq]l] with V(x) = — [sin(27rx) + .25 sin(47rx)]/(27r) as the ratchet potential. We focus on 
the region of parameters where, in |l8| (for q = 1), separation of masses was found. We fix 
7 = 2,KT = .1, r = .75, and D = .1875 and consider the values of the masses m = mi = 0.5 



and m = ni2 = 1-5 as in [18]. Our main result here is that the separation of masses is enhanced 



when a non-Gaussian noise with q > 1 is considered. In Fig. 4. a. we show J as function of q 
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for mi = 0.5 and m,2 = 1.5. It can be seen that there is an optimum value of q that maximizes 
the difference of currents. In that figure, this value, which is close to q = 1.25, is indicated with 
a vertical double arrow. Another double arrow indicates the separation of masses occurring for 
q = 1 (Gaussian OU forcing). The calculations are for a load force F = 0.25. We have observed 
that, when the value of the load force is varied, the difference between the curves remain ap- 
proximately constant but both are shifted together to positive or negative values (depending on 
the sign of the variation of the loading). By controlling this parameter it is possible to achieve, 
for example, the situation shown in Fig. 4.b., where, for the value of q at which the difference 
of currents is maximal, the heavy "species" remains static on average (has J = 0), while the 
light one has J > 0. Also it is possible to get the situation shown in Fig. 4.c, at which the two 
species moves in the opposite direction with equal absolute velocity. 

Summarizing, we have systematically studied the effect of a colored non Gaussian noise source 
on the transport properties of a Brownian motor. What we have found is that a departure from 
Gaussian behavior, given by a value of q larger than 1, induces a remarkable increase of the 
current together with an enhancement of the motor efficiency The latter shows, in addition, an 
optimum value for a given degree of non Gaussianity. When inertia is taken into account we 
also find a considerable increment in the mass separation capability. 

It is worth mentioning that in the studies of the influence of non Gaussian noises in stochastic 
resonance and other related phenomena, the system's response enhancement occurs for q < 1, 
as discussed in JTT| , 12, p!3| . The reason can be traced back to the dependence of the mean-first- 
passage-time with q, and the interplay between transition rates and modulation frequency. In 
contrast, in the present work, and due to the "favorable" influence of high values of the noise for 
increasing the transport effect, we observe that the relevant effect occurs in the opposite case, 
that is for q > 1. 

These studies could be of interest for possible relation with biologically motivated problems 

§ 



10, 21] as well as for potential technological applications, for instance in " nanomechanics" 



More specific studies of aspects in one or other area will be the subject of further work. 
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Figure 1: Current (a) and efficiency (b) as functions of q. The solid line corresponds to the 
analytical results in the adiabatic approximation while the line with squares shows results from 
simulations. All calculations are for m = 0, 7 = 1, KT = 0.5, F = 0.1, D = 1 and r = 100/(27r). 



Figure 2: Current (a) and efficiency (b) as functions of q. Results from simulations at KT = 
for D = 1 (circles), D = 10 (squares), and D = 20 (triangles). All calculations are for m = 
0, 7 = 1, F = 0.1 and r = 100/(2tt). 



Figure 3: Current as function of the mass for different values of q: circles q = 1 (Gaussian noise 
case), squares q = 1.3. 



Figure 4: Separation of masses: results from simulations for the current as a function of q for 
particles of masses m = 0.5 (hollow circles) and m = 1.5 (solid squares). Calculations for three 
different values of the load force: F = .025 (a), F = 0.02 (b) and F = 0.03 (c). 
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